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We derive an effective spin Hamiltonian for the one-dimensional half-filled asymmetric ionic Hub¬ 
bard model with alternating on-site interaction in the limit of strong repulsion. It is shown that 
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nating next-nearest-neighbor and three-spin couplings in the presence of a uniform and a staggered 
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I. INTRODUCTION 


During the last decades the correlation-induced metal-insulator (Mott) transition has been one of the challenging 
problems in condensed matter physics^ In most cases the translational symmetry is broken in the Mott insulator^ 
A notable exception is the one-dimensional repulsive (U > 0) Hubbard models 

'H'Hub t ^ ^ ^ ^ (^nQ^n+l,Q "J~ ^n+l,a^na) ^ ^ ^ Pn^Pn^ (^) 

n 0 = 7 , 4 , n 

at half-filling, where the dynamical generation of a charge gap is not accompanied by the breaking of a discrete 
symmetry^ In equation Q we have used standard notation, namely c^ a ( c na ) for electron creation (annihilation) 
operators and p na = c' na c na for the particle density at site n with spin projection a. The exact solution of the model 
CD in the case of a half-filled band reveals that the ground state is uniform, with exponentially decaying density 
correlations^ At the same time, spin excitations are gapless and thus magnetic correlations decay only algebraically^ 
This is readily understood in the large -U limit: indeed, for U |f | the infrared behavior of the model (Hj) at half-filling 
is fully described by the S'I/(2)-symmetric spin 5 = 1/2 Heisenberg Hamiltonian 


n Hei s = J^s n -s„ +1 + j'J2 s n-s 


n+2 5 


( 2 ) 


where J = 4 1 2 /U — 16 t 4 /C/ 3 and J' = 4 t 4 /U 3 up to the fourth-order terms in t/U Z& Since the condition |f| <C U 
implies that the frustration is weak J 1 -C J, the next-nearest exchange is irrelevant and the low-energy behavior of 
the initial electron system is governed by the standard isotropic Heisenberg Hamiltonian TLHeis = J^2 n ^n ' S n +i- 
Elegant techniques have been developed for calculating higher-order corrections to the Hamiltonian @. These terms 
are also irrelevant and leave the featureless character of the ground state intact 

The spin sector may even remain translationally invariant in the case of an explicitly broken translational invariance 
of the electronic Hamiltonian. For example, let us consider a scenario where two types of atoms are located respectively 
on even and odd sites of the lattice, with different on-site energies and/or different on-site couplings for the electrons. 
The Hamiltonian of such an extended version of the Hubbard model is given by 

n = t J2 ( c ™ c »+l,« + C i+l, a C na ) + UJ2 ( 1 + (-ir^PntPnl + y X Pna , (3) 

n,a n n a 


where 0 < <5, A/U <C 1. It possesses spin SU( 2) symmetry, but the translational symmetry has been reduced due to 
the doubling of the unit cell. At <5 = 0 and A/0 this Hamiltonian corresponds to the ionic Hubbard model (IHM)^ 
where electrons on even and odd sites have different on-site energies ±A/2, while at A = 0 and d / 0 equation © 
represents the alternating -U Hubbard model; 13 , where the electrons experience different on-site interactions on even 
and odd sites. 

At U = 0 the half-filled ionic Hubbard model describes a regular band insulator with equal charge and spin gaps 
and a long-range ordered (LRO) charge-density-wave (CDW) in the ground state. With increasing U the system 
undergoes two phase transitions, a first one at U = U c i from the CDW-insulator to a LRO dimerized insulator, and 


2 


a second one at U = U C 2 > U c \ from the dimerized phase to a strongly correlated (Mott) insulator— At U = U c 2 
the spin gap vanishes and the low-energy behavior of the system for U > U c 2 is again described by the Heisenberg 
Hamiltonian ED, with the difference that the spin exchange parameters J and J' now weakly depend on A. The 
broken translational symmetry of the model manifests itself only in the charge degrees of freedom via the presence 
of a LRO CDW pattern which persists even in the limit of strong repulsion, with the amplitude approaching zero at 

U -*■ 00 H 

The weak-coupling renormalization group analysis of the repulsive alternating -U Hubbard model (A = 0 and 
17(1 ±5) > 0) shows a qualitatively similar low-energy behavior at half-filling as the usual Hubbard model. Scattering 
processes arising from the alternating part of the interaction, which are relevant in the commensurates case of 1/4- and 
3/4-filled bands^ are irrelevant at 1/2-filliiig where the properties of the system are governed by the uniform part of 
the interaction. In the limit of strong on-site repulsion (U |t|), the infrared behavior of the alternating -U Hubbard 
model is once again described by the Heisenberg Hamiltonian ((2D , but with a slight modification - an alternating next- 
nearest-neighbor (nnn) exchange )T) r) [j/ + (—l) ra J[] S ra • S, l+2 — Numerical and analytical studies of the Heisenberg 
chain with alternating nnn exchange show that in the pertinent case of weak frustration (J 1 -C J), the alternation 
of J' is irrelevant and the infrared behavior of the model is fully described by the standard Heisenberg model with 
nearest-neighbor exchange ^ 7 Thus, even though the Hamiltonian 0 describes a fermion system on a lattice with 
broken translational symmetry, the information about the unit cell doubling at half-filling is fully accommodated 
within the high-energy degrees of freedom; the low-energy behavior of the system is described by a translationally 
invariant, isotropic spin Hamiltonian. It has to be noted that the above conclusion does not remain valid in the 
presence of bond alternation, i.e. if the hopping amplitude t is replaced by to + (—In this case one obtains that 
in the strong-coupling limit at half-filling the effective Hamiltonian, still given by the Heisenberg model, contains an 
alternating nearest-neighbor exchange ]T) n [Jo + (—l)”7i] S„ • S n +i, which leads to the spin-Peierls instability with 
gapped spin excitations— 

We now turn our attention to a model having full translational symmetry, but explicitly broken spin SU ( 2 ) sym¬ 
metry, the so-called spin-asymmetric Hubbard Hamiltonian 

^ ^ ^ 7" ^n+l 7 a^na) 7" ^ ^ , Pn*\Pn , (4) 

n,ot n 

where the hopping is spin-dependent (tf ^ t±). This model, introduced in the early 1990si* to interpolate between 
the standard Hubbard model (t- j- = fj,) and the Falicov-Kimball models ( t -j- > 0,tj, = 0), has been intensively studied 
during the last two decades— r— Away from half-filling the spin-up and spin-down particles are segregated in the 
ground state for large enough repulsion, both for the Falicov-Kimball model 27 and for the Hamiltonian ED with 
tf / tj, ^ 0— Therefore the spin-asymmetric Hubbard model appears to be well suited for studying transitions 
between phase-separated and homogeneous states, especially in one dimension— r— 

More recently, increased interest in low-dimensional correlated fermion models with spin-dependent hopping has 
been triggered by the fascinating progress in experimental studies of low-dimensional mixtures of optically trapped 
ultracold atoms of two different types,— such as ultracold atoms loaded into spin-dependent optical lattice s 33 ' 34 or 
trapped atoms of different masses , 35 ' 36 The great freedom available for generating optical lattices has also allowed 
one to play with the lattice geometry and to create bipartite lattices, which turned out to be a key ingredient 
for achieving higher-band condensatescoherence control^ density-wave dynamics^ and even graphene-like 
physics , 42 ' 43 It has to be emphasized that mixtures of fermions with different hopping amplitudes naturally appear in 
solid-state systems as well, namely when several bands cross the Fermi surface. This happens for instance in mixed- 
valence materials, organic superconductors^ small radius nanotubes^ and even graphene-based heterostructures— 
However, experiments with trapped ultracold atoms can actually engineer quantum many-body states and thus realize 
models of correlated fermions and bosons which are not available in usual solid-state structures— Recent theoretical 
predictions of various unconventional superfluid or superconducting*^— insulatin g 55 ' 56 and magnetic^ phases in 
such novel systems have further stimulated the interest in the spin-asymmetric Hubbard model. 

The broken SU{2) spin symmetry of the model (4) at / fj, is manifestly seen for a half-filled band in the 
strong-coupling limit (U |t-f|, \ti\), where to leading order the infrared behavior of the system is described by the 
anisotropic XXZ Heisenberg Hamiltonian 

n = J E ( S n S n +1 + S n S n+l + 7 ^ +1 ) , (5) 

n 

with J = 4t^ti/U and 7 = (f 2 + As the anisotropy parameter |^y| is larger than 1 for arbitrary t- j- 7 ^ tj,, 

the system has a finite spin gap and long-range antiferromagnetic order in the ground state— Nevertheless, the 
translational invariance of the initial lattice model ED is retained by the effective Hamiltonian ©, even if the ground 
state has lower symmetry due to the general phenomenon of spontaneous symmetry breaking. 
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In a recent paper we have studied the one-dimensional spin-asymmetric ionic Hubbard model in the limit of 
strong on-site repulsion (for a half-filled band)^ We have shown that for tf ^ the doubling of the unit cell by 
the alternating ionic potential A ^ 0 directly manifests itself in the spin degrees of freedom, and the effective spin 
Hamiltonian in the strong-coupling limit is given by the anisotropic XXZ Heisenberg chain with a staggered magnetic 
field 


n = jJ2( S nSn+l + ^^1+7%) -hJ2(-D n S Z n 


where 


J = 


4 Ut 




U( 1-x 2 ) ’ 


7 = 


tt + tj 
2 t t t t 


h = 


2 (t 2 — t 2 )a; 

U(l-x 2 ) 


( 6 ) 

(7) 


and x = A/U. For tf ^ ij, and finite x, the translational symmetry is broken already at the level of the effective spin 
Hamiltonian via the presence of the staggered magnetic field. Since this represents a strongly relevant perturbation 
to the spin system, the ground state is characterized by a long-range antiferromagnetic order with explicitly broken 
translational symmetry. The excitation spectrum is gapped and the gap exhibits power-law dependence on the 
parameter h. 60 

In the present paper we extend our analysis to the case of explicitly broken translational symmetry in the on-site 
interaction and derive the effective spin Hamiltonian for the one-dimensional spin-asymmetric alternating -U ionic 
Hubbard chain represented by 

K = J2 ta ( C L c n+l, a + 4+l,a C na) + U ( 1 + ME' 5 ) P„fPnl + J Pna ■ (8) 

n,a n n a 


We find that up to fourth-order terms in t a /U , the infrared behavior of the lattice fermion model © at half-filling is 
governed by the following effective spin Hamiltonian: 


Keff = E J±(S%SZ +1 + sy n s y n+1 ) + J\\(S*S* +1 - £) + J' ± (n) (S^ +2 + Sy n S y n+2 ) + J'^n) (S* n S* +2 - i) 

n 

+ E \w±(n) [(S^s- + ! + SlS y n+1 )\ + W\\ (n) - E h{n) S* , 

n L -I n 

where 


J± = 

J\\ = 

Ji(n) = 
J||(n) = 

W_i_(n) = 
W\\ (n) = 




1/(1 — A 2 ) 
2 (*? + *?) 


1 - 


2 (tf 




U 2 (l - A 2 ) 2 
6 (tf + tf) 


2A 2 - 


U( 1-A 2 ) U 3 (l — A 2 ) 3 


(1 + 3A 2 ) - 


1-A 2 
1 — <5 2 

4 tpf 


U 3 ( 1 - A 2 ) 3 


5-A 2 - 


4(1 - A 2 ) 
1-5 2 


4 tftf 


17 3 (1 -A 2 ) 3 

2 (tf + tf) 

U 3 (l - A 2 ) 3 

+ (-l) 


2 + 2A 2 — ^ A !/ l +(-!)" 


1 + 3A 2 + 
25 


1 — 5 2 
1-A 4 
1 — A 2 




4 tpl 


t/ 3 (l — A 2 )(l - A 2 ) 
1 + 5 2 


1/3(1 -A 2 ) 2 1 — <5 2 


1/ 3 (1 — A 2 ) 2 (l — 5 2 ) 


[4 tftf — (tf + tf)( 1 + A 2 )], 


4Atftf(t 2 — t 2 ) f 25 
t/ 3 (l-A 2 ) 2 \ 1 — A 2 

4A (tf-tf) ( 25 


+ (-l) 


+ 


1 — <5 2 


1/3(1-A 2 ) 2 (1-d 2 


+ (-l) r 


3 +A 2 
l^A 2 
3 +A 2 2 
1 - A 2 + 1 - 6 2 


h(n) = h 0 + (-l) n h\ = 

+ (-l) 


2A 5(tf-tf) 


1/3(1 - A 2 ) 2 (l — 5 2 ) 
2 A (tf - tf) 


+ 


1/(1-A 2 ) 


1 - 


t 2 4 - 1 2 

+ ^ 


2t/ 2 (l-A 2 ) 2 


5(3 +A 2 ) - 


2(1 - A 2 ) 1 
1 -5 2 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 
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with A = 5 + A/U. 

As we observe, the leading terms (oc I/ -1 ) are the same as in Q except that the parameter x is replaced by A. The 
higher-order terms (oc I/~ 3 ) include the renormalization of the nearest-neighbor coupling J, the next-nearest-neighbor 
exchange with an alternating part whose existence is determined solely by 5, and corrections to the magnetic field 
alongside the less conventional three-spin terms, all having both homogeneous and alternating parts. We also note 
that the expressions of the fourth-order terms obtained in our earlier worb^S are not entirely correct and they should 
be replaced by the appropriate limit (5 = 0) of the above results. 

A detailed derivation of the expressions (!§!)- (1TT>1) is presented in the following. In Sections HT1 and EH a unitary 
transformation is applied to the electronic Hamiltonian in the case of a half-filled band, eliminating hopping processes 
between many-electron states with different numbers of doubly occupied sites. In Section IIVI we briefly discuss the 
Hubbard operators, which are used in the subsequent Section El to derive the effective spin Hamiltonian. Finally, 
Section fVTl summarizes the main results of the paper, while the |Appendix| contains some technical calculations of the 
spin exchange terms. 


II. THE STRONG-COUPLING APPROACH 

In the strong coupling limit (U |f|), the perturbative treatment of the half-filled Hubbard model based on 
expansion of the Hamiltonian in powers of t/U goes back to the original derivation of the effective spin Hamiltonian 
to the second order by Anderson^i Afterwards, using different versions of the degenerate perturbation theory, effective 
spin Hamiltonians up to higher orders in t/U have been obtained. In particular, Klein and Seita§£ derived the sixth- 
order spin interaction for the Hubbard chain, while Bulaevskii and Takahashi^ obtained the fourth-order terms for 
the half-filled Hubbard model in higher dimensions. More recently, these perturbative methods have also been applied 
to Hubbard models with more general interactions . 63 

An alternative approach to construct the effective Hamiltonian is based on unitary transformations. Harris and 
Lange^ used such a transformation to obtain second-order results and to calculate spectral properties of the Hubbard 
model. A transformation which systematically incorporates higher orders in t/U has been proposed by Chao, Spalek, 
and Oles^ In their expansion, closed expressions for the effective spin exchange are obtained to any order. However, 
beyond the second order their method is not very well controlled since the transformation of the Hamiltonian involves 
an approximation for the band energies, and higher-order terms mixing different Hubbard bands are not eliminated 
properly . 66,67 

A consistent transformation scheme which allows one to remove all unphysical terms and to derive the t/ [/-expansion 
up to any desired order has been formulated by MacDonald, Girvin and YoshiokajS In their scheme, interaction terms 
which do not conserve the number of local electron pairs are eliminated from the Hamiltonian order by order in an 
iterative treatment, generating new interactions and thus improving the accuracy of the transformation at each step. 
Later their approach has been successfully employed to obtain effective spin Hamiltonians in the case of extended 
versions of the Hubbard model on a square lattice with next-nearest- and next-next-nearest-neighbor hoppings . 68,69 

Another consistent scheme for construction of the effective spin Hamiltonian up to any given order in powers of 
t/U has been developed by Stein^ who utilized Wegner’s method 70 of continuous unitary transformations with 
subsequent solution of the corresponding flow equations for the half-filled Hubbard model. Later a similar approach 
has been used to reveal an additional (hidden) symmetry of the Hubbard model on any bipartite lattice ^ 1 

In this paper we apply the method developed by MacDonald, Girvin and Yoshioka for the standard Hubbard model*! 
to the Hamiltonian T-L = T + V, where 

T = taC n<* C ma » ( 17 ) 

<n,m> a 

V = ~2 A-ta + U° P 2 ra-H.tA 2 n-H.-L + Pln.fPln.i. > (1®) 

n,ot n n 

and the brackets in the sum <n,m> signify that n and m are labels for neighboring sites. The on-site couplings 
U 0 = U{ 1 — S ) and U e = U( 1 + 5) are supposed to be strong, U e > U 0 |£-f|, \t±\, A, implying that the parameters 
5 = ( U e — U 0 )/{U e + U 0 ) and A = 8 + A/U satisfy the conditions 0 < <5 < 1, 0 < A < 1. 

In the larg e-U limit of the standard Hubbard model (5 = A = 0) the many-electron states are grouped according 
to the number of doubly occupied sites (doublons) Nj- In the present case with S, A > 0 these Hubbard subbands 
are split into groups of states classified by two numbers, N,i e and N,i 0 , representing the numbers of doubly occupied 
sites on even and odd sublattices, respectively. The hopping operator T mixes the states of these subbands. The 
“unmixing” can be achieved by introducing suitable linear combinations of the uncorrelated basis states. The S 
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matrix for this transformation, and the transformed Hamiltonian, 

W = e lS Ue- tS , (19) 

are generated by an iterative procedure, which results in an expansion in powers of the hopping amplitudes fy and/or 
ti divided by the on-site energies U e and/or U 0 - 

This expansion is based on the separation of the kinetic part of the Hamiltonian into three terms: 

T = T 0 + Ti + T_i, (20) 

where To leaves the number of doubly occupied sites unchanged, and T\ (T_i) increases (decreases) this number by 
one. In the present case of broken translational symmetry each of these contributions is further split into several 
different terms, depending on whether the electron hops from an even to an odd site or vice versa. 

In particular, the To term is split into four separate processes: 

T 0 = Tq e + Tg ° + T<f + T<f . (21) 


Here 


and 


ir= E E*°a 

<2n,m> a 


P 2 n,o) ^2 n,a^ma 


(! - Pma ) 


( 22 ) 


^0 — X] X] ta (1 P2n+l,a) C 2n+l,a C ma (1 Pmoi ) (23) 

<2n+l,ra> a 

correspond respectively to hopping processes where an electron with spin a hops from a singly occupied odd (even) 
site to an empty neighboring even (odd) site, while 

X) T.^A^maPm, ( 24 ) 

<2n,m> a 


and 


rpdo _ 

-*0 — 


z^ z^ 

<2n+l,m> a 


P2n-\-l,a ^2n+l,a^ma Prt 


(25) 


represent hopping processes where an electron with spin a hops from a doubly occupied odd (even) site to a neighboring 
even (odd) site which is already occupied by another electron with the opposite spin a. 

In a similar fashion, the operators T± i, which change the number of doublons by one, are also separated into even 
and odd parts T ±1 = TJ/ + T/ 1; where 

^ ^ y tg p2n.a^2n.a^ma Pma) (26) 

<2n,m> o; 

and 

■^1 ^ ^ ^2n+l,a ^2n+l,a^ma Pma) 

<2n+l,m> a 

increase the number of doublons on the sublattice of even (odd) sites, while 

1 ^ y ^ ^ Pnd) ^na^2m,a ^2m,a (^) 

<n,2m> a 

and 

r ^'— 1 ^ ^ (1 Pno) ^na^2m+l,o: ^2m+l,a ’ (^) 

<n,2m+l> o: 

respectively decrease the number of doublons on the even and odd sublattices. 
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One can easily check the following commutation relations: 

[v, t‘] = (p + <y Ml 0 )A a r*, 


where fi = 0, ±1 and 


A s = 


A, 

-A, 

(U e - U 0 ) + A, 
~(U e - U 0 ) - A, 
U e +A, 
U D - a, 


s = pe 
s = po 
s = de 
s = do 
s = e 
s = o 


(30) 


(31) 


The relations (l3()l) reflect the fact that the energy of the system changes by (/x + 0 )A S as a result of the hopping 
process T®. 


III. EFFECTIVE HAMILTONIAN IN THE CASE OF A HALF-FILLED BAND 


Let us now search for the unitary transformation S which eliminates hops between states with different numbers 
of doubly occupied sites in the transformed Hamiltonian 

W = e iS Ue~ iS =n+ [iS,H] + ^[iS, [iS,H]} +.... (32) 

We follow a recursive scheme^ which allows to determine such a transformation to any desired order in t a /U. The 
last two terms of the initial Hamiltonian 


H = H ,(1) = V + T 0 + Tr + T_i 


(33) 


may be transformed away by choosing 


iS = iS w 


1 (T?-T°_ 1 ) + -^—(T?-Tl 1 ). 


U 0 - A 


U e 


Substituting and m into the expansion (E2D and applying pi)) . we obtain 


H '{2) =V + Tq + 

+ 

+ 

+ 


1 


U 0 - A 
1 


U e T A 

U e +Uo 


{[T^Toj + lTo^ + iT^T^]} 
{[T^Tol + lTcTfJ + ^T^]} 

{ [T® ,T° 1 ] + [T °, T®-J } 


2(H e + A)(Z7 0 - A) 
U 0 -U e - 2A 
2(E7 e + A)(L/ 0 - A) 


{ [T®, T°] + [TfT® J } + 0(t 3 /U 2 ). 


(34) 


(35) 


We focus on the case of a half-filled band, where in the large-C/ limit the lowest-energy states I'Ll) have exactly one 
electron at each site. In this subspace no hops are possible without increasing the number of doubly occupied sites. 
Therefore, 


T! 1 |'Ll}=0, T° 1 |^l)=0, T o |4-l)=0, 


and the effective Hamiltonian (1551) is reduced to 


n' {2) 


T° X T? 
U 0 - A 


TVP 

U e + A 


o(t 3 /u 2 ). 


(36) 


(37) 


To proceed further, we define: 




( 38 ) 
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Using we can write 


V,T^[{s}, { M }]1 = E \ Si (m + ^,o)T (l) [W, {^}] ■ 


(39) 


i =1 


77'^^ contains terms of order (t a ) k , denoted by 'H’' lk \ which couple states in different subspaces. By definition 
[U, ^ 0 and can be expressed in the following way: 


K ,[k] = E E c {a } (M)r (fc) [{«}, M] 

{a} { M } 


^ ' Mi 7^ 9 . 


If at each fc-th order step we choose <S^ a s 
where 5^1 is the solution of the equation 


g( k ) = + s 

[isw,v] = - n ,[k] 


and therefore equals 


(I ,,ti 

A w ’ r<*)[{<■},M], 

{ Q } { m } Z^i=i A ai(.M* + <Vi,o) 


^ ' Mi 7 ^ 0 . 


i=l 


then the transformed Hamiltonian 




/(fc+i) _ e iS^ ne -iS^ 


(40) 

(41) 

(42) 

(43) 

(44) 


contains terms up to the order of (t a ) k /U k ~ 1 which commute with the unperturbed Hamiltonian V and mix states 
within each subspace only. 

The conditions (15C>1) can be generalized to higher orders 


if 


r (fe) [W,{M}] |*l> = 0 , 


^p[M] = E ^ <0 

i=p 


(45) 


(46) 


for at least one value of p. Equation (14511 can be used to eliminate many terms from the expansion for T-L' in the 
subspace of minimal (V). 

The final expression of the transformed Hamiltonian W up to the fourth order reads: 

rpo r nO rpe rpe rpo rppOrpperpo rpo rpdorpderpo rpe rpde'Y'do'J^e 

o,/(4) _ 1 -l 1 1 _ _-1^1 _ 1 - 1 1 0 1 0 1 1 _ 1 - 1 1 C) 1 0 1 1 _ 1 0 1 l 

U a - A U e + A U 0 (U 0 - A) 2 U e (U 0 — A) 2 U a (U e + A) 2 

pP^ppOpe po pP°pdope pe pdepP e pO 

~ U e {u e + A) 2 " {Uo - A)U a (U e + A) “ (t/ 0 - A)t/ 0 (U e + A) 

po rpdorpPO rpe pe pP e pdepO po pe rporpe 


(U Q - A)U e (U e + A) (Uo- A)U e (U e + A) (U e + A){U 0 - A)(U e + U a ) 


pe rpo rperpo 


pe rpo rporpe 


p^o pe pepo 


(U e + A)(U 0 -A){U e + U 0 ) (U e + A) 2 (U e + U 0 ) (Uo- Af(U e + U 0 ) 

po po popo pe pe pepe po popo po pe pepe pe 


2 (Uo - A)3 2 (Ue + A) 3 (Uo - A)3 ( U e + Af 

U 0 + U e 


2(U 0 -A) 2 (U e + Ay 


[ r o iT o T e iT e + T e_ lT e T ° lT °] . 


(47) 
























IV. HUBBARD OPERATORS 


To handle the effects of strong interaction properly, it is important to know whether at the beginning or at the end 
of a given hopping process a particular site is doubly occupied or not. For this purpose one introduces the so-called 
Hubbard operators^ X“ b , which are defined at each site of the lattice and describe all possible transitions between 
the local basis states |a), | b): unoccupied |0), singly occupied with an up-spin It) or a down-spin |t) electron, and 
doubly occupied |2). The original electron creation (annihilation) operators can be expressed in terms of Hubbard 
operators in the following way: 


c t = X 


0:0 


rj{a) X, 


2ol 


Cna = X"“ +rj(a) X , 


a2 


(48) 


where rj(a) = { ^ = J[ 

Conversely, in terms of creation (annihilation) operators the Hubbard operators have the form: 

X a0 = c 1 (1-b-) 

na\ x ricxJ » 

Y““ = c f c - 

n '"not na ’ 

X 0 n 0 = (l-p nt )(l-p nl ), 

X n a = PnM ~ Pna) ■ 

The Hubbard operators containing an even (odd) number of electron creation/annihilation operators are Bose-like 
(Fermi-like) operators. They obey the following on-site multiplication rule 


X n a = y(a) 4aPia > 

X^® = 71 (ol) C f -ct 
yv n l l\ Ui ) '"naFna j 

Y 22 _ - 

Pn\ Pn\. ’ 


(49) 


X™X™ = 5 qr X* s 


(50) 


and commutation relations 


[X”X%\ ± = 5 n 


± 6 ps x %), 


(51) 


where the upper sign stands for the case when both operators are Fermi-like, otherwise the lower sign should be 
adopted. 

It is straightforward to represent the hopping terms introduced in Section QI] by the Hubbard operators: 


Tg° 

= E 

V'' 4- Y a 0 yC )cx 

2.^ L a yv 2n+l yv m > 

TZ e = 


< 2 n+l,m> cx. 


'Y'do 

1 0 

= E 

V' + Y^cx Y a 2 

L cx yv 2n-\-l yv m > 

II 


< 2 n+l,ra> cx 



= E 

E ri(a)t a XiZ +1 X° m ", 

II 


< 2 n+l,ra> 

CX 


T°_ x 

= E 

£ 77 (a) t a X«°X|2 +1 , 

,1 

II 


<n,2m+l> a 


<2n,m> cx 

- E E 

<2n,m> ot 

E E^ 

<2n,m> ol 

- E E 

<n,2m> cx 


r 2oc yot.2 
^2n yv m ? 


(52) 


One also easily verifies that the X-operators describing the transitions between singly occupied states can be rewritten 
in terms of spin 5=1/2 operators as 


Vl~l = J r = ,Q+ Y4T = J r — Q- 

n u n > '-'njAn'|~ ‘-’n i 

Xtt = I + 5~, X^=i-5^. 


(53) 


V. THE SPIN HAMILTONIAN 

Using the relations (I5l21) -(l53 l) . it is straightforward to rewrite the products of T-terms in (HTl) via the Hubbard and 
hence the spin 5 = 1/2 operators. We first consider the simplest two-component T-terms at great length to elucidate 
the procedure for more complicated contributions. 


A. The second-order terms 

Let us start from the hopping term which corresponds to creation and subsequent annihilation of a single doublon 
on an even 2n-th site. Since the electron hopping is restricted to nearest neighbor sites, this process only includes 
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electrons located on two neighboring sites 2n ± 1, and is given by: 


T+T-f 


E \ \ A j.2 y ctO y c*2 y2a y Oct 

/ v v c a ^ 2 n+q^- 2 n^ 2 n^- 2 n +9 


+ +- yo:0 ya2 y2a: yOa: 
L a L a y ^2n+q^2n^ L 2n^2n+q 


L/2 


E V ' \ ^ <2 vaa vaa / / yaa yaa 

/ v / v c a ^-271+9^-271 c a+a ^2n+q uX 2n 


n—1 q—zLl a 
L/2 


j .2 ytt y - 14 - I 4.2 y -14 ytt _ 

6 t ^-271+9^-271 ' U ^-271+9^-271 

1 \ / 1 \ . 0/1 


Ml ( 


ytt ytt I ytt ytt 
^ 271 + 9^-271 ■ ^ 271 + 9^271 


E E ( 

71=1 9=±1 
L/2 

EE[*t(J + Si n+q ) (i - S|„) +t| (± - Sl n+q ) (± + Sf„) 

('^’ti+g+n + +n++ 2 ri) 

E E [ (4 + *1) (J - 5|„5| n+9 ) - \ (*? - tl) (Sf„ - Sf„ + J 

-2i t +Sf„Sf„ +9 + 5| n 5|„ +g )l 


71=1 9=±1 

L/2 

E 

71=1 9=±1 


E [ - 2 * 7 * 4 . (S*S * +1 + S%S% +1 ) - (4 + tj) (S*S * +1 - 1) - (- 1 )" (4 4) + 


(54) 


As to the second term TZ\T°, which describes creation and annihilation of a pair on an odd site, the calculation is 
essentially the same, with the only difference being that one needs to make a replacement 2n f> 2n + 1. Below this 
shift is absorbed in q: 

L/2 

rpo rpo \ A \ ^ \ ^ j.2 yaO ya2 y 2a yOa J. j. _ yaO y<52 y 2a yOa 

1 -l 1 1 — 2_^ 2_^ 2_^ I c a ^ 271 + 9 ^ 271 + 1 ^- 271 + 1 ^ 271+9 L a L a ^2n+q^2n+l^2n+l^2n+q 
7i=l 9=0,2 a 


E E [4(l + s L +q )(l-SL + x) + 4(l- s L +q )(\ + SL + i) 

7i=l 9=0,2 

— Ml {^2n+q^2n+l + +ra++2n+l) 

L/ 2 . 1 

EEl^^Hr SL +g SL + 0 + 2 (4~tl) (SL +q SL+4 

n=l 9 = 0,2 

—2t^ti(S 2n+q S 2n+ i + 

^2n+q^2n+l)\ 

L 

E [ - 2 * t +S*S* +1 + S£S£ +1 ) - (4 +t|) (S*S* +1 - -) + (-1)” (4 - $ S* 

71=1 


(55) 


The combination of these two processes yields the second-order effective spin Hamiltonian: 


where 


' 1 /( 2 ) _ _ 

f f — 


L eff 


U„ — A 


t o_ iT o _ 


1 


U e + A 


T+T-f = 


</E^ 


71+1 


■ « +1 + 7 (S*S* +1 - 1/4)] - hJ2(-l) n S z n , 


4<t*4. + * 1 . 2 +t - ffl 

U(1 — A 2 ) ’ 7 2t t t4. ’ C/(l — A 2 ) ' 


(56) 


(57) 


As we see, the second-order effective Hamiltonian, which describes the spin degrees of freedom of the initial lattice 
fermion model, is the Hamiltonian of spin S' = 1/2 frustrated XXZ Heisenberg chain in the presence of a staggered 
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magnetic field. The amplitude of this field is proportional to the product of the parameter A quantifying the broken 
translational symmetry of the underlying fermion model, and the spin-dependent hopping asymmetry parameter 
tf — ip Thus, in contrast with the spin-isotropic case (t- 1 - = t±), the infrared properties of the spin-asymmetric model 
are described by a Hamiltonian with an explicitly broken translational symmetry. 

It is instructive to check several limiting cases. In the case of spin-symmetric electron hopping (t-j- = = t), the 

effective Hamiltonian (15H1) reduces to the Hamiltonian of the isotropic (5'C/(2)-invariant) Heisenberg chain 

n eff = J^Sn -Sn+1, (58) 

n 

with a uniform exchange constant J = At 2 /U(l — A 2 ). Thus, even if the translational symmetry of the underlying 
fermion model is broken (A / 0), the second-order effective spin Hamiltonian remains translationally invariant. 

In the complementary case of the Hubbard model with spin-dependent hopping (A = 0, iq ^ fq), the second-order 
effective Hamiltonian properly reflects the broken spin symmetry and is given by the Hamiltonian of anisotropic 
(H(l)-invariant) Heisenberg chain 


H eff = JJ2 ( S n S n+l + « +1 + J S^ +1 ) 


(59) 


with the anistropy parameter 7 = (t 2 + t^)/2t^ti > 1^2 

Finally, in the limiting case of the Falicov-Kimball model (t± = 0), the second-order effective spin Hamiltonian 
reduces to the Ising model in a staggered magnetic field: 


^e 2 // = E( J i \s z n s* n+1 -(-ir h s~), 

n 


where 


2 f 2 2 At 2 

J|1 = U(l - A 2 ) ’ h= U(1 — A 2 ) ' 


(60) 


(61) 


The physical mechanism responsible for appearance of the staggered magnetic field in the effective spin Hamiltonian 
(1561) can easily be traced in the ultimate limit of the Falicov-Kimball model, however the argument remains valid 
also for arbitrary f-j- > t± > 0. Due to the doubling of the lattice unit cell, energetically it is preferable to locate all 
immobile fermions on odd sites, while the mobile up-spin fermions will predominantly occupy even sites. In this limit, 
the process of creation and annihilation of a doublon takes place only on odd sites and gives rise to the following 
Ising-type spin exchange parameter Jy 1 ^ = 2t 2 /(U 0 — A), while in the opposite case, where all immobile spins are 

located on even sites, the same process yields the exchange constant = 2 t 2 /(U e + A). The difference between the 
exchange energies for these two patterns equals 




4Af 2 

U{ 1-A 2 ) 


= 2 h. 


(62) 


B. The fourth-order terms 

The same technique as the one employed in the previous section can be used to rewrite the products of four T-terms 
in the effective Hamiltonian (ITT!) via the spin S = 1/2 operators. There are 18 terms of this type. It is convenient to 
unite them in groups characterized by the similarity of the hopping processes and by the number of created doublons 
at the intermediate steps. 


1. Group A: Four-T product terms of the form T-\TqTqT\ 

There are eight terms of this type in the effective Hamiltonian (1471) . In these processes the number of created doubly 
occupied sites is one. Four terms correspond to processes where the doublon is created and eventually annihilated on 
the same site, while the other four terms describe processes where the doublon is created on an odd (even) site and 
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annihilated on the neighboring even (odd) site. The calculations are straightforward and one obtains the following 
expressions for the operators (the details can be found in the | Appendix!) : 


f 4 — f 4 


T— 1 Tq°Tq & T° = £ [ - tMA + 4)( S Z S Z+i + S V nS V n+i) - 2 t 2 t tl(S z n S* n+1 1/4) , g 

n 

~~ (4 — tfy {^2n-1^2n+l ~ 1/4) ~ 2(t| — ^) — 1/4)S'|„ +1 — 


(-1)"S* - 


-2^(4 - i|)[(Sf n _iS£, + S|„_iS!JS| n+1 + 5|„_ 1 (5f„5f„ +1 + S'|„5 2 V 1 )] 


(63) 


Tl.TfTfT? = E [ - *t*4(*t + *!)(^+i + «+1) - 2i?iJ(S*S* +1 - 1/4) + ^±(-1)"^ - 

n 

— (*t — ^4-) (‘^2n‘S'2n+2 _ 1/4) + 2(fc|. — t±) ('Sfra^fra+l _ 1/4) + 

+ 2*T*t(*T - <l)[(^^+l + «„ +1 )Sf n +2 + (^2n+l ^2n+2 + ^n + l^ +2 )] 

T^Tl e Tl°Tl = E[-*t*l(*t + + S^+i) - 2^^(5';S'; +1 - 1/4) - ^±(-l)"S£ - 

n 

~~(4 — ^) 2 (‘^2n‘S'2n+2 - 1/4) ~ 2(f 4 — t|) ('S'fn'^fn+l - 1/4)5271+2 — 

— 2*t*4-(*t — ^4.) [(^271*^27^+1 + ‘S , 2n5'2ri+l)‘5 , 2n+2 + (^ra+l ^'2n+2 + Sfra+l £2/1+2)] ! 

T^TtTtn = E[-*4*4$ + ^)(5*S*+1 + «+i) - 2^f (£*£*+1 - 1/4) - ^=^(-1 )”£* - 

n 

-(t? - i|) 2 (^n-i^+i - 1/4) + 2(4 - tj)(£f n _i££„ - 1/4)£f„ +1 + 

+2hh(# - tU^n-iSL + SL- 1 SL)SL +1 + SL-i(s% n s$ n+ i + s y 2 n s y 2n+1 )] 
Tt x Ttirn = E [ - *r*4(*?+*!)( s n^+i+«+i) - 2^ (s*s* +1 -1/4) + 

n 

+itfti (4 ~ i|)[ {^ 2 n- 1 ^ 2 n ~ & 2 n- 1 ^ 2 n) ~ (‘S’fn'^ln+l — ^ra^fra+l) ] + 

+ 21^1/ (§2ro—1 ' §2ra+l — 1/4) , 

T— r Tq C Tq b T° = E [ - *4*4$ + + ^n^n+i) - 244 (£*£* +1 - 1/4) - 

n 

—it^(4 ~ i|)[ {^2nS 2 n +1 ~ ^n^n+l) — (£fn+l £271+2 — ^2n+1^2n+2) \ + 

+ 2fcftJ (S2n ' S2 n +2 — 1/4) , 

T— 1 Tq°Tq°T^ = E [ - *r*4(*t + fyffiSZ+i + «+i) - 2 44 (S*S* +1 - 1/4) + 

n 

+it^tl(4 ~ i|)[ {S 2n ^2n+1 ~ ^n^n+l) ~ (£?ra+l $271+2 ~ ^n+l^n+2) ] + 

+^44 (§2ra ' S2„+2 — 1/4) , 

TZiTq°Tq°Ti = E [ - + t D( S n S n +1 + 5 ^n+l) - 2i?i| (£*£^x - 1/4) - 

n 

—it\tl(4 ~ ^f)[ {^2n-1^2n ~ ^2n-1^2n) ~ (‘^2ra‘S'2n+l ~ ^fn^fn+l) ] + 

+ 21fl/ (S2n-1'82™+! — 1/4) 


(64) 


(65) 


( 66 ) 


(67) 


( 68 ) 


(69) 


(70) 


2. Group Bl: Four-T product terms of the form TlxT^LiT^Tf and T!Li T“iTiT“ (a 7 ^ &/ 

There are four terms of this type in the effective Hamiltonian (S3. These terms correspond to processes where 
two doublons are created on even and odd sites and then consecutively annihilated. If the first pair is created on site 
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£, the hopping processes leading to the creation of the second pair are restricted by the existence of an empty site 
adjacent to £. Consequently, using the relations (ITT^ll — (foTTll one obtains a restricted double summation of the form 


T 0 _ x T e _ x T°Tl = 

n,m^n,n±l 


- 2 t^(s* s* +1 + sy m si +1 ) (t? + tl){S* m S* m+1 - 1/4) 


K-l ) m (S^-S^ +1 ) 


- 2 + SlS y n+1 ) - (tf + - 1/4) - 

= (TZ x Tf) * (T® 1 Tf). 


£tzC ( _i)™ ( ^-^ +i ) 


(71) 


Here we have introduced the notation * to denote multiplication of infinite sums over the indices n and m with the 
restrictive condition m ^ n, n ± 1 . 

For the other terms of the same group we analogously obtain 


T e iT O iT o T e = (TO iT o) * (T^T?) , 

ri 1 r° 1 T 1 e T 1 ° = tz x tZiT x t° = (Tf ^f) * (r^rf) . 


(72) 

(73) 


However, since all the hopping processes in these products take place on disjoint pairs of sites (n, n + 1) and (to, to + 1), 
one can freely commute the iS-operators past each other, so that the order of the multiplicands becomes irrelevant: 


(TZiT°) * [TZ X TZ) = (TfiTf) * (TZiT?) = - [(TZ X T?) * (Tf X T^) + (T^Tf) * (TZiT°)\ . 


(74) 


3. Group B2: Four-T product terms of the form T'“ 1 T“ 1 T“Tf 


There are two terms of this type in the effective Hamiltonian (1T71) . These terms correspond to processes where two 
doubly occupied sites are created either on even or on odd sites and then consecutively annihilated. As before, the 
creation of the first pair puts limitations on the processes responsible for the creation of the second pair. In addition, 
since there are two different ways how one can get the same configuration corresponding to the pair of doublons 
located on two odd or two even sites, an extra factor of 2 appears in the expressions for these terms: 


rj-UD rj-iO rjnOrj^O _ 

= 2 E 

n,m^n,n±l 


- 2t t h( S m S m+l + £'m<S'm+l) “ (*t + i |)0 S m' S m+l 

- 2 t t t±(szsz +1 + sy n s y n+1 ) - + tl)(s*s * n+ 1 - 


-i/ 4 ) + ^(-ir(^-^ +1 ) 

1/4) + %^(—!)"(£'- S* +1 ) 


= 2(Tf jTf) * (T^Tf), 


(75) 


TZxTZiTfTf = 


= 2 E 

n,m^n,n± 1 


- 2tMS x m S x m+ 1 + S V mS V m+1 ) - m + tl)(S* m S* m+1 - 1/4) 

- 2 i t i±(S*S * +1 + S y S y +1 ) - C 1 1 + tJ)(S*5 * +1 - 1/4) - 


_ _^1/ 1 \m(az _ qz \ 

2 V L ) \°m D m+l) 

^(-i nsz-sz+j 


= 2(TZ 1 TZ) * (Tf 1 T 1 e ). 


(76) 


4- Group C: Four-T product terms of the form T_iTiT_iTi 

There are four terms of this type in the effective Hamiltonian (TI71) . These terms correspond to processes where a 
doublon is created and immediately annihilated on a site £ and then another doublon is created and annihilated on 
an arbitrary site £'. Using ©-(El), we obtain 

To_ iT o T o iT o = E _ 2 tMS x m S x m+l + sy m s y m+1 ) - + tl)(S* m S z m+1 - 1/4) + ^(-l) m (^ - S* m+1 ) 

n,m 

- 2 Mt^S*S* n+1 + S y S y +1 ) - C 1 1 + tl)(S z n SZ +1 - 1/4) + 3=£(-l)"(S* - S~ + i) 

= {TZ X T?) ■ (TZ X T °), 


(77) 
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'pe pepe pe 


E 

n,m 


2t th( S m S m+l + ^-S'l+l) - (i| + ^)( S m S m + l - 

• ' - 2t x t i {S*S* +l + SlSl +1 ) - (t 2 + tf)(S* n S* n+1 


l/4)-^_i!(-l ) m (S* m -S* m+1 ) 
-l/4)-^(-l TO-5^) 


= (T^Tf) ■ {T^Tf) 


( 78 ) 


and 


Tf 1 T 1 °Ti 1 T 1 e = (Tf ^f) • (TfiTf), 


T e iT e r o iT o = ( T e iT e) . _ 


(79) 


As we observe, in marked contrast with the first eight terms ©-(EH) which only couple spins located on neighboring 
sites, the remaining ten terms given by CD-© contain countless number of all possible two-spin S^S^, three-spin 
5 'n 5 n+i 5 'm and four-spin SPS p n+1 S^S q m+1 combinations, where p,q = x,y,z. The situation is rescued by the fact 
that after combining identical terms in the Hamiltonian (|47|) . each term of the type (TfjTf) * (T^ 1 Tf) will have its 
counterpart of the type (T“ X T“) • (T h _ l T l l ‘) with just the opposite coefficient. As a result, all terms corresponding to 
distant spin-spin interaction are canceled: 


= E 

n,m=n,n± 1 


( T « iT o) . (J’° 1 yo) _ ( T o iT o) * = 

- 2 tMS x m S x m+ i + S^^+i) - (*? + - 1/4) + %^(-l) m (^ - 5^ +1 ) 


• [ - 2 tMs x n s x n+l + s q s v n+1 ) - (t 2 + t 2 )(s*sz +1 -1/4) + Vk-i T(SZ - S z n+i ) 
E { [ - 4W*t + *?)(^ +1 + « +1 ) - (3 4 + 2 t 2 t 2 + 3tJ)(5*5* +1 - 1/4)+ 


+2t?t2(S*5 * +2 + S»S^ 2 ) + (*| + <t )(^,^ +2 - 1/4) + (-If 


+4_+4 

t c 4- 


(45*_ 1 5*^ 1 +55*) 


+2t t *4(*t - W-iOM+i + ^^+i) + + S^SM+i] ] } , 


(80) 


and 


E 

n,m=n,n=tl 


( T e iT e) . ( T e iT e) _ (ye ^e) * (Tf ^f) = 

- 2^(5- S * +1 + S* S£ + i) - (t 2 + t 2 )(S^ m+1 1/4) - 3z£(-l)"*(S* - ^+i) 

• [ - 2 tMS x n S x n+ i + S« n S y n+ 1 ) - (if + tl){S* n S z n+1 - 1/4) - ^(-1)"(S* - S* +1 ) 
E{[-4*t*4(*t + ^X^+i + « +1 ) - (3if + 2 ifif + 3$(S*S * +1 - 1/4)+ 

n ^ L 


+ 2 ifif(S*S* +2 + sy n s y n+2 ) + (if + tJ)(S*S* +2 -1/4) - (-ir 


-(4 S’^S’S; 


n+1 


■5S*) 


+2*r*4(*t - *!)[^-i( 5 n^ + i + S"SJ[ +1 ) + + ^_i^)^ + i]] } 


(81) 


2 [( T -i T i°) • (TlrTf) - (Tf 1 T°) * (Tf^f)] + - [(T® X T®) • (Tf ^ - (Tf^f) * (Tf iTf)] = 


E { - + tl)(S x n S x n+ 1 + ^^ +1 ) - 8ifif (S*S* +1 - J) + 

n 

+2t 2 t 2 (S* n S* n+2 + S"S" +2 ) + 2tftJ(S*S * +2 - i)} . 


(82) 


For the last two terms we have taken their combination to avoid calculation of extra terms which cancel each other 
in the sum. 

Inserting finally the relations (IMll - dH^ll into (IT71) we obtain that up to the fourth order the strong-coupling effective 
spin Hamiltonian for the spin-asymmetric alternating-U ionic Hubbard model is given by 


rpo rpc 

0 /( 4 ) _ _ 

H eff ~ 


Tf 1 T 1 e 


U 0 - A U e + A 
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where 


1 

' rpo rppOrpperpo rpderpdorpe " 

1 -1 1 0 1 0 1 1 . 1 -1 1 0 1 0 -*1 

1 

" rpo rpdo'T'de'T'O rpe rpperppOrpe~ 

2 -1-^0 -*0 1 1 , ^-l^O 1 0 1 1 

4o 

L Mo -A) 2 (4e + A) 2 J 

" 4e 

[ (Uo- A) 2 (C/e + A) 2 J 


1 


(Uo - A )U 0 {U e + A) 

1 


| 'rj-iO rppOrj-ydorpe _|_ rj~ie rj~\derj-iperj-\0 j _ 

+ T e iT pe T de T oj + 


(t/ 0 - A)U e {XJ e + A) 

_ 1 _ [ (T O iT o } . (T^Tf) _ (T O iT o) * (T^T”)] + 

' [(TfiT-f) • (T^Tf) - (TfiTf) * (T^Tf)] + 
1 


(Ue + A) 3 

C/o + C/ ( 


(U a - A) 2 (U e + A) : 


-{ 2 [( T -i T i°) ' M-iM) - M-iM) * M-iM)] + 


1 


+ - [(T^T-f) • (T° r T°) - (T^Tf) * (TfiTf 


E K^+i + « +1 ) + -J)]-E MM ^ 

n n 

+ E [ J i(M ( S n S n+2 + S V n S y n+2 ) + J(j(n) (S*S * +2 - i)' 

n 

+ E Kw [(Sn-iSn + + S 2 _i(« + i + « +1 )] + W„(«) 


(83) 


J_l — 


Uo- A Ue + A 


+ MlMf + 4))' 


+ 


— 4tfti(t^ + t 2 ) 


U 0 U e l(U Q - A) 2 (U e + A) 2 (U 0 - A)(U e + A) 

1 1 M 


L([/ 0 - A) 3 (14 +A) 3 (U 0 — A) 2 (U e + A) 2 


t/ 3 (1 — 6 2 ) (1 - A 2 ) 2 


J ii — 


i , hhi 1 t + t V 2 

U 1 - A 2 

4 Ml 

C/(l — A 2 ) 

(*?+$ 


EE*? " " 2 


1 - 


2 (4 + 4) 

4 2 (1 — A 2 ) 2 

1 1 


2 + 2\ 2 - 


1 - A 2 


1 -<5 2 


+ 2 ^ 2 ^° + ^ 


Uo - A U e + A 

1 


-3(4 + 4) 

1 


1 


U 3 (1 - A 2 ) 3 


1 


(Uo- A) 3 (t/e+A) 3 J 

2 


4^ L([/ 0 -A) 2 (t/ e + A) 2 (U 0 — A)(U e + A)J 


- 2 tit 


2+2 


t 4 

2,4-2 


+ 


4 (U 0 + U e ) 


(Uo- A) 3 (U e + A) 3 (J7„- A) 2 (4e + A) 2 


2(£ 2 + i 2 ) 6(i| + i|) 


MM = (- 


17(1 - A 2 ) ?7 3 (1 - A 2 ) 3 

1 


(1 + 3A 2 )^ 


+44 


t/ 3 (l-A 2 ) 3 


!)"{(* 


2 -M 


4 4 4[/ o _ a [/ e + aJ 


+ 


(4-4) 


i 


5 — A 2 — 


1 


4(1 - A 2 ) 
1 - <5 2 


2 (U e + Uo) - (-l) n (U e - Uo) l(Uo - A) 2 (Ue+ A) 2 

+4 j.4 


EiEv r_ 1 


i 


= (-i) ? 


2 l(U 0 - A) 3 (U e + A) 3 
J 2A(1 2 — l 2 ) A(4 — tf) 


U( 1 - A 2 ) ?7 3 (1 - A 2 ) 3 


5 (3 +A 2 ) - 


2\ 2(1-A 2 ) 


1 — S 2 


2 A 5(4 - 4 


(84) 


(85) 


^(l - A 2 ) 2 (l - <5 2 ) ’ 


( 86 ) 
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J±(n) — 2A/i 


T 4 l(U e + Uo) + (-1 )"(E/ e - U 0 )](U 0 - A )(U e + A) 

1 1 U e + Uo 


2 t 2 A 


+ 


+ 


V(U 0 - A) 3 (U e + A) 3 (Uo- A) 2 (E/ e + A) 2 J 




C/ 3 (1 - A 2 ) 3 


2 + 2A — 


2 ( 1 -A 2 ) 


2x21 


1 — A 2 


+ (-l) r 


4 ^t*l 


J,|(n) = 


(+2 _ . 2x2 _“_ 

(E/ e + Uo) + (~l) n (U e - Uo) 
4 


U 3 (l - A 2 )(l — 5 2 ) ’ 

1 1 


- 2tlt 2 


'-(Uo — A) 2 (C/e + A) 2 J 
+ 


T A [(E/e + E/„) + (-l)"(E/ e - C/ 0 )](C/ 0 - A)(C/ e + A) 


+ (*t + *t) 


1 


1 


L(E/ 0 -A)3 (U e + A) 3 J 


2 /^ 


c/ e + c/ D 


4^(1 + < 5 2 ) 


+ 


2(/f + 


1 + 3A 2 + 


(Uo - A) 2 (U e + A) 2 
1 — A 4 " 1 


1 — <5 2 


C/ 3 (l - A 2 ) 2 (l - S' 2 ) U 3 (1-A 2 ) 3 
+ ( — -^"rrsTT-12X271-Hx [4/f/J - (/f + i l)( 1 + A 2 )] , 


+ 


C/ 3 (l — A 2 ) 2 (l — 5 2 ) 
W±(n ) = (-l)"{2t t ^(i|-t2 ) 


1 


1 


+ 2.t^t±(t 2 — t 2 ) 


(U e + Uo) - (-1 ) n (U e - Uo) V (Uo - A) 2 (E/e + A) 2 

1 1 
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( 88 ) 


(89) 


(90) 


It should be pointed out that the presented expressions (1551) - (TO1) fully agree with the known results in the limiting 
cases of the standard^ (tf = t± = t, A = 5 = A = 0) and the alternating-C/— (tf = ty = t, A = 0, A = 5 ^ 0) 
Hubbard models, but some of the fourth-order coefficients do not coincide with the results obtained previously for 
the ionioi^ (t^ = = t. 5 = 0, A = A/C/ = x) and the spin-asymmetrio^ (tf / tj., A = 5 = A = 0) Hubbard 

models. More specifically, for the ionic Hubbard chain we arrive at a different expression in the numerator of the 
nearest-neighbor coupling J, whereas for the spin-asymmetric Hubbard model the disparities concern the numerators 
of the coefficients Jy and . We presume that these discrepancies are due to the perturbative schemes adopted by 
the authors of Refs. — and — , as it is known that some of these procedures are not sufficiently well-controlled at 
higher ordersj 66 i 67 


VI. CONCLUSION 

We have derived the effective spin Hamiltonian for the low-energy sector of the one-dimensional half-filled spin- 
asymmetric alternating -U ionic Hubbard model in the limit of strong on-site repulsion. The obtained Hamiltonian is 
that of a frustrated Heisenberg chain with alternating next-nearest-neighbor exchange and three-spin coupling in the 
presence of a uniform and a staggered magnetic field. As expected, the nnn exchange is larger for two spins separated 
by a site with low on-site repulsion than for spins separated by a site with high on-site repulsion. The intensity of 
the three-spin coupling and the amplitudes of the magnetic fields are proportional to the product of the parameter 
A, which reflects the broken translational symmetry of the lattice, and the difference between up- and down-spin 
electron hopping amplitudes f-|- — tp The most dominant effect however comes from the staggered magnetic field, 
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and therefore, in marked contrast with the spin-isotropic case t- j- = t±, the ground-state properties of the considered 
electron system are described by a spin-chain model with explicitly broken translational symmetry. 

We also remark that the general picture outlined above remains valid in the case of a half-filled bipartite lattice of 
a higher dimension - to the lowest order one again obtains the anisotropic nearest-neighbor spin exchange and the 
staggered magnetic field which, as before, dominates any higher-order terms arising from the more complex lattice 
geometry. 
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Appendix: Derivation of the effective spin exchange expressions for the four-T product terms of the form 


T-iTqTqTi 


In order to rewrite the products of four T-operators in the language of spin S = 1/2 operators, once again we make 
use of the Hubbard X-operators defined in section El the procedure is in essence the same as the one employed in 
section IV Al for the terms consisting of two T-operators. 

Let us consider for example the term T^T^T^T^, which belongs to the group of the processes where the electron 
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pair is created and annihilated on the same site. In the summations below, the square brackets around the lattice 
indices [n, to, k\ indicate that for an even site 2 to its odd partners 2k + 1 and 2n +1 represent neighboring sites. Thus, 
for a fixed to we have two possible sets: k = m — 1 and n = to, or k = m and n = m — 1 . 
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(91) 


In a similar manner one can derive the expressions of the operators T^iT^T^Tf, T^T+T+Tf and Tf y+T+T'f, 
belonging to the same group as the operator considered above. However, it should be noted that one can obtain all 
of these expressions directly from (1911) by switching the roles of the odd and the even sites and/or interchanging the 
hopping amplitudes tf ++ : 
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Tl 1 Tg e Tg°T? = T^T^T^T? (t t o tj.) = 
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The expressions for the remaining four operators, corresponding to the processes where the pair is created on one 
site and annihilated on a neighboring one, can be established in an analogous way (utilizing where necessary the 
freedom of renaming the lattice indices k -O- n to facilitate the calculation): 
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